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A  SENSING  ARRAY  SYSTEM  WITH 
IMAGE  STATISTICS  PROCESSING 


INTRODUCTION 

A  semiautomated  technique  for  extracting  a  selected  set  of  cartographic  features 
such  as  road  intersections,  straight-line  roads,  and  rectangular  objects  from 
aerial  photographic  imagery  using  Walsh  transform  was  recently  reported.1,2 
An  effort  to  extract  and  recognize  another  set  of  selected  area  cartographic 
features  such  as  water,  forest,  city,  field,  and  farmland  with  a  spatial  domain 
approach  is  described.  The  purpose  of  this  work  was  to  investigate  and  verify 
the  performance  of  an  automated  apparatus  for  feature  extraction  using  available 
off-the-shelf  components.  This  apparatus  can  be  used  for  features  from  photo¬ 
graphic  as  well  as  radar  imagery.  The  images  which  were  used  in  this  study  were 
aerial  photographs  taken  from  two  locations. 


Among  many  available  image  spatial  domain  processing  techniques,  the  image 
histogram  and  the  image  texture,  which  make  use  of  first  and  second  order  image 
statistics,  respectively,  were  considered  for  implementation.  The  hardware  system 
consists  of  a  32-  by  32-element  solid-state  array  to  convert  cartographic  images 
into  electronic  signals,  a  minicomputer  to  process  the  electronic  signals  from  the 
array  using  a  controlled  software,  and  a  computer-controlled  two-dimensional 
translation  stage  as  the  image  holder.  The  image  statistical  processing  technique 
is  implemented  as  a  software  for  the  system. 


A  brief  system  description  is  followed  by  a  discussion  on  the  selection  of  a 
feature  vector  from  both  image  histogram  and  image  texture  methods.  Although 
many  components  of  the  feature  vector  were  calculated  for  each  incoming  test 
image,  only  four  components  were  used  for  classification  purposes.  Classification 
was  based  on  two  types  of  classifiers  with  optical  (aerial)  imagery  as  test  input. 
It  was  found  that  both  classification  methods  resulted  in  a  recognition  accuracy 
of  better  than  90  percent.  Finally,  conclusions  are  given. 


lP.r.  Chen  and  W.W.  SeeniuUer,  Signal  Signature  of  Topographic  Features  Using  Analog  Technology. 
ETL  Report  No.  0185,  U.S.  Army  Engineer  Topographic  Laboratories,  I  ort  Belvoir,  Va.,  May  1979. 
AD-A076  1 10. 

^J.R.  Singleton  and  P.L.  Chen,  Detecting  Line  Hoad  and  Roaddntersection  Patterns  at  Various  Angles. 
U.S.  Army  Engineer  Topographic  Laboratories,  1  ort  Belvoir,  Va.,  ETL-0274,  October  1981. 
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SYSTEM  DESCRIPTION 

The  block  diagram  of  the  system  is  shown  in  figure  1 .  A  9-  by  9-inch  trans¬ 
parency  is  illuminated  by  a  white  light  source,  and  a  section  of  the  image  is 
projected  onto  a  Reticon  32-  by  32-element  solid-state  array  through  an 
imaging  lens.  The  array  converts  the  optical  energy  of  the  image  into  a  video 
signal.  The  video  signal  is  quantized  into  10  bits  of  digital  signals,  and  the 
latter  are  sent  to  the  Hewlett-Packard  2108  minicomputer  for  processing.  The 
computer  first  takes  in  the  quantized  signals  of  32  by  32  pixels  of  1,024 
gray  levels  array.  With  the  brightest  and  darkest  pixels  within  a  frame  as  the 
maximum  and  minimum,  this  quantized  image  array  is  next  scaled  down  to 
become  32  by  32  pixels  of  16  gray  levels.  The  histogram  and  joint  prob¬ 
ability  matrix  are  then  obtained  using  the  scaled  array  as  the  input.  The  next 
step  is  to  compute  a  feature  vector  based  on  the  histogram  and  the  joint  prob¬ 
ability  matrix  separately.  Although  16  components  of  the  feature  vector  were 
computed,  only  4  components  are  used  for  classification  of  the  selected  carto¬ 
graphic  features  from  optical  imagery.  Either  a  sequential  template-matching 
classifier  or  a  minimum-distance  classifier  is  used  to  classify  the  input  images 
into  one  of  the  preassigned  image  categories,  or  it  is  recognized  as  a  reject  if  it 
does  not  belong  to  any  of  these  categories.  The  classification  result  is  then  in¬ 
dicated  on  the  CRT  console.  At  the  end  of  classification  a  signal  is  sent  to  the 
translational  stage  controllers  to  move  the  stages  in  the  predetermined  \  and  v 
positions,  and  a  new  section  of  image  is  projected  onto  the  surface  of  the  solid- 
state  array.  This  procedure  repeats  itself  until  all  preselected  image  sections  are 
classified.  The  software  for  the  classifiers  is  listed  in  the  appendix. 


FIGURE  1.  System  Block  Diagram. 
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FEATURE  VECTOR  SELECTION 


Many  publications  are  available  on  the  image  statistics  processing  techniques.* 
The  most  concise  form  of  the  feature  vector  based  on  the  first-  and  second- 
order  image  pixel  amplitude  distribution  was  given  by  Pratt.3  He  defined  a 
discrete  image  array,  F(j,  k),  and  the  first-order  probability  distribution  of 
image  amplitude  of  a  measurement  window  centered  about  (j,  k)  as 


P(b) 


N(b) 

M 


0) 


where  M  represents  the  total  number  of  pixels  in  the  measurement  window, 
N(b)  is  the  number  of  pixels  of  amplitude  b  in  the  window,  and  0  <  b  <  L  -  1 , 
and  L  is  the  number  of  gray  levels  of  FG,  k).  The  following  measures  have 
been  formulated  by  Pratt  as  a  concise  means  of  describing  the  shape  of  first- 
order  image  histograms:4 


L  -  i 

Mean  b  =  £  b  P(b)  (2) 

b  =  0 


L  -  i 

Variance  a2  =  £  (b  -  b)2  P(b)  (3) 

b  =  0 


Skewness  bs 


L  —  1 


E 

b  =  0 


(b  -  b)3  P(b) 


(4) 


Kurtosis  bk 


1  L-l 

—  E  (b  -  b)4  P(b)  -  3  (5) 

ab  b  =  0 


*See  Pratt,  Haralick,  Sutton,  and  Duda  in  the  Bibliography  for  sources  of  information  on  this  topic. 
3W.K.  Pratt,  Digital  Image  Processing,  John  Wiley  and  Sons,  Inc.,  New  York,  1978. 


L-l 

Entropy  bE  =  —  £  P(b)  log2  [  P(b)]  (?) 

t  =  o 


The  second-order  histogram  features  are  based  on  the  definition  of  the  joint 
probability  distribution  of  pairs  of  pixels.  Haralick  et  al.  have  proposed  a  number 
of  texture  measures  based  on  the  joint  amplitude  histogram  of  pairs  of  geometri¬ 
cally  related  image  points.5  Pratt  stated  that  the  two-dimensional  histogram  can 
be  considered  as  an  estimate  of  joint  probability  distribution.  Consider  a  pair 
of  pixels  F(j,  k)  and  F(m,  n)  that  are  separated  by  7  radial  units,  at  an  angle 
d  with  respect  to  the  x-axis  of  the  measurement  window.  The  histogram  estimate 
of  the  second-order  distribution  is  given  by  Pratt  as6 7 


P(a,  b) 


N(a,b) 

M 


(8) 


where  M  is  the  total  number  of  all  occurrences  in  the  measurement  window,  and 
where  N(a,  b)  is  the  number  of  occurrences  for  which  F(j*  k)  =  a.  Fun,  n)  =  b. 
The  following  are  the  texture  measures  that  were  used  in  this  study:7  8 


L-i  l  —  1 

Mean:  a  =  £  £  a  P(a,  b) 

a  =  0  b  =  0 


L  -  1  L  —  1 

b  =  £  E  bp(a’b^  (io> 

a  =  0  b  =  0 


5R.M.  HaiaJick,  k.  Shanmugan,  and  I.  Dinstein,  Texture  Features  for  Image  Classification,  II  I!  Trans¬ 
actions  on  Systems,  Man,  and  Cybernetics.  SMC-3.  November  1973. 

^W.K.  Pratt,  Digital  Image  Processing.  John  Wiley  and  Sons,  lnc„  New  York.  1978. 

7Ibid. 

®RJd.  HaiaUck.K.  Shanmugan,  and  I.  Dinstein.  Texture  Features  Jor  lma.ee  Classification.  11  I  1  Irans- 
actions  on  Systems,  Man,  and  Cybernetics,  SMC-3,  November  1973. 
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(11) 


Variance:  Va  =  £  E  (a  -  a)2  P(a,  b) 

a=  0  b=  o 


L  —  1  L  —  1 

vb  =  E  E  (b-b)2  P(a,  b) 

a=  0  b  =  0 


(12) 


L  —  1  L  —  1 

Covariance:  C0  =  £  £  (a  -  a)  (b  -  b)  P(a,  b)  (13) 

a=0  a=0 


L  —  1  L  —  1 

Autocorrelation:  Au  =  a  b  P(a,  b) 


a  =  0  b  =  0 


L  —  1  L-l 

Absolute  Value:  Ab  =  E  E  I  a  —  b  I  P(a,  b) 

a=0  b=0 


L-l  L-l 

Energy  :  Eg  =  £  £  [P(a,  b>l: 

a  =  0  b  =  0 


L" 1  L“L  P(a,  b) 

Inverse  Difference:  Id  =  V  V]  - 

a  =  0  b  =  o  1  +  (a  -  b)2 


L-l  L-l 


Inertia:  In  =  £  £  (a  -  b)2  P(a,  b) 

a  =  0  b  =  0 


(14) 


(15) 


(16) 


(17) 


(18) 


L-l  L-l 

Entropy:  En  =  -£  £  P(a-  b)  log2  (P(a,  b))  (19) 


! 


f 
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For  our  application  the  joint  probability  matrix  was  made  to  be  symmetrical 
so  that  a  =  b,  and  Va  =  Vb.  The  variation  of  P(a,  b)  due  to  different  values 
of  6  was  not  treated  in  this  report,  since  0  =  0  degree  provided  sufficiently 
separable  feature  vector  for  all  the  image  categories  considered.  Each  input 
image  was  first  scaled  down  to  1,024  pixels  of  16  gray  levels  (L  =  16).  Fifteen 
components  for  the  feature  vector  based  on  equations  (1)  through  (19)  were 
then  computed.  It  was  discovered  that  only  a  measure  from  the  first-order  histo¬ 
gram,  two  measures  from  the  second-order  statistics  (texture),  and  a  measure 
of  the  number  of  array  elements  above  the  threshold  value  (NAEATV)  were 
needed  for  classification  purposes. 


CLASSIFIER  DESIGN 

For  the  selected  set  of  imagery,  only  three  components  of  the  feature  vector 
computed  in  the  previous  section  plus  an  NAEATV  measure  were  required  to 
constitute  a  4-dimensional  classifier.  These  components  of  the  feature  vector 
are  as  follows: 


1.  Skewness  (Histogram) 

2.  Covariance  (texture) 

3.  Autocorrelation  (texture) 

4.  The  number  of  array  elements  above  the 
threshold  value  (NAEATV) 


Two  types  of  classifier,  a  sequential  template-matching  classifier  and  a  minimum- 
distance  classifier,  were  chosen  for  this  purpose. 


Sequential  Template-Matching  Classifier.  In  order  to  determine  the  upper 
and  lower  limits  of  the  template  values  for  each  image  category,  many  pro¬ 
totypes  were  obtained  from  imagery.  Figure  1  chows  the  template  ranges  for 
each  component  of  the  feature  vector  selected  for  classifying  the  selected  set 
of  imagery.  Six  template  ranges  for  the  covariance  were  designated  as  follows: 
water,  -1.0  to  2.25;  field,  0  to  3.5;  forest,  1.0  to  4.0;  city,  4.0  to  9.0;  farm¬ 
land,  9.0  to  20.0;  and  not  recognized,  above  20.0.  The  template  ranges  for  the 
autocorrelation  were  assigned:  forest,  0  to  90;  field,  70  to  300;  and  water, 
90  to  300.  For  discriminating  “forest”  from  “field,”  the  template  ranges 
of  the  skewness  were  chosen  as  field,  -10  to  0,  and  forest,  above  0  to  5.  Like¬ 
wise,  the  template  ranges  for  the  NAHATV  were  set  to  be  water,  0  to  300 
and  field,  301  to  1023.  It  is  noticed  that  many  template  ranges  overlap  one 
another:  however,  with  the  sequential  template-matching  scheme  discussed  next, 
most  of  the  images  tested  so  far  have  been  recognized  or  classified  correctly. 


The  computed  covariance  of  the  unknown  incoming  input  image  is  first  compared 
to  the  template  range  of  the  covariance.  If  it  is  within  the  range  of -1.0  to  2.25, 
it  is  temporarily  assumed  to  be  in  the  category  of  “water.”  The  next  step  is  to 
examine  whether  the  autocorrelation  of  the  unknown  image  is  in  the  range  of  90 
to  300  or  not.  If  it  is  in  that  range,  the  NALATV  of  the  unknown  image  is 
tested  against  its  corresponding  template  range.  If  it  falls  in  between  the  range  of 
0  to  300,  the  unknown  image  is  classified  as  “water.” 


If  any  one  of  the  previously  described  tests  is  negative,  the  unknown  image  is 
temporarily  assumed  to  be  in  the  category  of  “field.”  Just  as  in  the  tests  per¬ 
formed  for  the  category  of  “water,”  the  covariance,  autocorrelation,  skewness, 
and  NALATV  of  the  unknown  image  are  now  sequentially  compared  to  their 
corresponding  template  ranges  (see  FDCN2,  appendix  A).  If  all  tests  are  true, 
the  unknown  image  is  classified  as  “field.”  If  any  one  of  these  four  steps  is 
false,  the  unknown  image  is  temporarily  reassigned  to  the  category  of  “forest.” 


Likewise,  the  unknown  image  has  to  go  through  three  sequential  tests  in  this 
assumed  category  for  its  covariance,  autocorrelation,  and  skewness.  The  strategy 
is  just  as  stated  before.  If  all  the  tests  are  true,  the  unknown  is  classified  as 
“forest”  or  it  is  reassigned  to  the  next  category  of  image,  and  so  on.  It  is  noticed 
that  only  a  single  covariance  comparison  is  needed  for  the  categories  of  “city” 
and  “farmland”  (see  FDCN2,  appendix  A).  Finally,  if  the  unknown  image 
does  not  belong  to  any  stage  of  the  tests  described  above,  it  is  then  classified  as 
“not  recognized.” 
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FIGURE  2.  Template  Ranges  for  Feature  Vector  Used  for  Classification. 


Minimum-Distance  Classifier.  Assuming  that  the  image  categories  of  interest 
(or  the  classes  to  which  the  input  imagery  are  to  be  classified)  are  represented 
by  prototypes  (or  reference  patterns)  Zl ,  Z2, .  .  . ,  Zm  ,  the  Euclidean  distance 
between  an  arbitrary  vector  of  unknown  classification,  X,  and  the  i-th  pro¬ 
type  is  given  by 

Dj  =  HX-Zjll  =  v/(X-Z;)T  (X-Z;)  (20) 


where  (X  -  Zj)T  is  the  transpose  of  (X  -  Zx).  A  minimum  distance  classifier 
computes  the  distance  from  a  pattern  X  of  unknown  classification  to  the  pro¬ 
totype  of  each  category  and  assigns  the  pattern  to  the  category  to  which  it  is 
closest.  In  other  words,  X  is  assigned  to  category  “  i  ”  if  <  Dj  for  all 

j  ^i* 


To  establish  the  prototypes  Z] ,  Z2,  .  .  .  ,  Zm  ,  the  components  of  the  feature 
vector  — “covariance”  and  “autocorrelation”  (computed  from  the  joint  pro¬ 
bability  matrix),  the  component  “skewness”  (evaluated  from  the  first-order 
histogram),  and  “NAEATV”  were  used.  They  are  expressed  as 

=  [  (COV)j, (AU)j,  (NAEATV);,  (SKEW);  ]T  (21) 


The  required  five  prototypes,  one  for  each  image  category  of  interest,  i.e.,  city, 
field,  farmland,  forest,  and  water,  were  obtained  by  using  the  mean  value  of 
these  four  components  computed  from  many  image  samples.  These  Zj  were 
then  used  to  compute  D,  and  classify  the  incoming  unknown  pattern  X* 


The  feature  vector  components  “autocorrelation”  and  “NAEATV”  are  approx¬ 
imately  10  and  100  times  larger  than  the  other  two  components  for  all  five 
image  categories  considered.  In  order  for  the  smaller  components  not  to  be  over¬ 
powered  by  these  two  large  components  in  classification  computation,  “auto¬ 
correlation”  and  “NAEATV”  were  divided  by  factors  of  10  and  100  respec¬ 
tively,  Likewise  the  corresponding  components  of  the  feature  vector  of  the 
unknown  incoming  pattern  X  are  also  normalized  in  the  same  way  before  being 
submitted  to  the  classifier  for  computation.  Figure  3  shows  the  final  com¬ 
ponents  for  the  prototypes  used  for  this  classifier. 
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FIGURE  3.  Prototypes  for  Minimum-Distance  Classifier. 


TEST  RESULTS 


A  set  of  high  quality,  1  to  20,000  scale,  aerial  photoimagery  from  Baltimore, 
Maryland,  and  Pittsburgh,  Pennsylvania,  consisting  of  image  categories  such  as 
city  (combination  of  commercial  and  residential  structures,  DLMS  category 
number  504  F1C  301  and  DLMS  category  number  505  F1C  401),  Field 
(agriculture  field  used  primarily  for  crop  and  pasture  lands,  DLMS  category 
number  510  F1C  950),  farmland  (cropland,  DLMS  category  number  510  FIC 
950),  forest  (deciduous,  DLMS  category  number  510  FIC  952),  and  water 
(river  and  inlet  to  bay,  smooth  fresh  water,  DLMS  category  number  501  FIC 
941)  was  used  for  this  experimentation.  Figures  4(a)  to  8(a)  show  the  line 
printer  output  of  the  typical  image  categories:  city,  field,  farmland,  forest,  and 
water,  respectively.  Each  is  printed  in  16  gray  shades  by  line-printing. 


In  Figures  4(b)  through  8(b),  the  16  components  of  the  typical  feature  vector 
computed  by  using  equations  1  to  19  and  the  “NAEATV”  are  shown  for  each 
image  category.  It  was  found  that  only  the  components  “covariance”  and 
“autocorrelation,”  evaluated  on  the  basis  of  image  texture,  and  the  “skewness,” 
computed  from  image  histogram,  are  relatively  separable  for  all  the  image  cate¬ 
gories  considered.  Therefore,  these  components  together  with  the  “NAEATV” 
were  chosen  to  constitute  a  four-dimensional  classifier  for  both  the  sequential 
template-matching  and  the  minimum-distance  classification  schemes. 


Approximately  120  image  sections  (for  all  5  classes)  were  randomly  scanned 
and  were  provided  as  the  input  test  imagery  for  evaluating  the  classification 
accuracy  for  both  classification  methods.  The  results  are  illustrated  in  Figures 
9  and  10.  Approximately  92  percent  overall  classification  accuracy  was  ob¬ 
tained  for  the  sequential  template-matching  classifier.  The  overall  classification 
accuracy  for  the  minimum-distance  classifier  is  relatively  higher  and  in  the  range 
of  about  97  percent.  The  results  presented  above  are  valid  for  high  quality  aerial 
photoimages.  The  images  were  taken  from  the  same  aircraft  altitudes,  but  covered 
two  different  geological  areas.  Various  parameters  such  as  image  quality,  scale 
changes,  resolution,  season  and  geological  formation  were  not  addressed. 


The  same  experimental  setup  was  used  to  extract  and  classify  a  set  of  radar 
imagery  with  modified  classifiers  (different  feature  vector,  template  range,  and 
prototypes).  Approximately  90  percent  recognition  accuracy  was  obtained  for 
this  selected  set  of  radar  imagery.  This  effort  will  be  reported  in  another  research 
note. 
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FIGURE  4.  (a)  Pictorial  Print  of  Input  Image,  and  (b)  Feature  Vector  for  City. 
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FIGURE  5.  (a)  Pictorial  Print  of  Input  Image,  and  (b)  Feature  Vector  for  Field. 
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FIGURE  6.  (a)  Pictorial  Print  of  Input  Image,  and  (b)  Feature  Vector  for  Farmland. 


i’8i  «».«  SECOND  ORDER  STA 


FIGURE  8.  (a)  Pictorial  Print  of  Input  Image,  and  (b)  Feature  Vector  for  Water. 


CLASSIFIED  CATAGORV 


CLASSIFIED  CATAGORY 


FIGURE  10.  Classification  Results  for  Minimum-Distance  Classifier. 


CONCLUSIONS 


1.  Statistical  processing  of  photographic  images  in  spatial  domain  provides  a 
convenient  means  for  evaluating  image  texture  and  coarseness  of  area  features. 

2.  The  technique  is  most  applicable  for  extracting  and  classifying  if  the  search 
window  contains  only  a  single  category  of  cartographic  features  such  as  city, 
forest,  field,  water,  or  farmland. 

3.  Better  than  90  percent  classification  accuracy  was  obtained  for  both  the 
sequential  template-matching  classifier  and  the  minimum-distance  classifier  for 
a  set  of  approximately  1 20  selected  aerial  photographic  images. 

4.  With  modified  feature  vector,  template  ranges,  and  prototypes,  the  method 
can  also  be  applied  to  extract  and  classify  radar  images. 

5.  The  scheme  may  produce  ambiguous  results  when  the  search  window  contains 
more  than  a  single  category  of  cartographic  area  features. 
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APPENDIX  A.  Computer  Printouts  FDCN2  and  FDCN4 


1=00004  IS  OH  CR00009  USING  00003  BLKS  R=0000 
F  1 M4  ,  L 


C 

SUBROUTINE  FDCN2  --  REV  03/20/8  1 

C 

c 

C  SUBROUTINE  10  PERFORM  FEATURE  CLASSIFICATION  FOR  PROGRAM 

C 

C  TXTR? 


C 

c 

SUBROUTINE  FOCN?(lTS.COV,AUSkW,LUOT  > 

1  F (  C  0  V  11  -1  OR  COV  CT  2  25  )  GOTO  010 

I  F  <  A  U  LT  90  OR  AU  CT  300  1  GOTO  010 
1  F (  ITS  I  T  0  OR  ITS  Cl  3  0  0  >GG1G  010 
UR  I T  E(  LUO  I  . 80  ) 

CU  T  0  8  08 

810  1F<C0V  LI  0  OR  COV  Cl  3  50 >  GOTO  820 

1  F  <  A  U  LT  70  OR  AU  GT  30  0  KOTO  8  20 
1FCSKU  LT  -10  OR  SLU  GT  0>  GOTO  820 
lh  ns  LT  3  00  OR  US  GT  1  0  24)  GOTO  820 
UR  IT  £(  LUO  T , 01  ) 

CO  TO  080 

820  I  F  (  COV  L  I  1  00  OR  COV  CT  4  0)  GOTO  810 

l  f  (  A IJ  LT  0  OR  AU  GT  90)G0  TO  830 
I  F  <  S  K  U  EO  0  OR  SKU  GT  5)  GOTO  830 
URUF  <  L  UOT  •  82  ) 

GO  TO  800 

830  lf(COV  L  I  4  0  OR  COV  G1  9  0)  GOTO  035 

UR  I TE<  LUO T , 03  > 

GOTO  886 

035  I F ( COV  L T  9  0  OR  COV  GT  20)COTO  840 
WRI  T  E<  LUO  T ,  85  ) 

GOTO  886 

84  0  UR  I  T  E<  LUO  T ,  84  ) 

GOTO  888 

90  F  OR  H  A  T  <  IX, ’WATER"  > 

81  F  OR  M  A  T (  1 X , ’FIELD’  > 

82  FORMAT (  1 X . ’FOREST’  ) 

83  F  OR  N  AT  <  IX. ’CITY  ’  ) 

84  F ORMAT (  I  X ,  ’ TH I S  CARTOGRAPHIC  FEATURE  IS  NOT  SPECIFIED’) 

85  F ORMAT <  I  X . ’ FARML AND  ’  ) 

888  1  F  4  L  UO 1  EO  6 > U R I T E < LU 0 T  ,  8 8 0  ) 

800  FORMAT ("  1  ’  ) 

RF  T  UPN 

END 

END! 
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APPENDIX  A.  Con  tin  u  ted 


4FPCH4  1*00004  IS  OH  CROOOO?  USING  00005  OIKS  R=0000 
0  00  1  f  T  114 

0002  [sSES*=55s»=a;= ja=r=  SUBROUTINE  F  0  C  M  4  as*tj*i9«a»*»s3«ai*sssss»*et»s;* 

000  3  C- 

0004  C -  HIM1HUH  DISTANCE  CLASSIF  ICO? ION  SUBFROCRAN 

0005  C  = 

0006  C  =  WRITTEN  04  /28/8  1 

0007  C *  RECOPIED  07/22/81 

0  00  s  c  = 

0009  Cs*=*=s--*===*-*:»as  =  ==  =  =  =  5,=  ;,=  =  =  s:a=!S==*R“"'  =  *»*3*»«*»*»«**»»»“:*»**»*** 

OOJO  SUBROUTINE  F D C M4 < CA V . AU ,  1 1 S , S K W , l UOI > 

0011  D I  HE  NS  ION  COVR<  7  >. ACORRt  7  >, I TSRRt 7  >, S*WR<  7 >, Dl S( 7  ) 

0012  DATA  COVR/7  197,0  095,13  70,2  906,0  36/ 

0013  DATA  AC0RR/40  76,181  2,90  03.65  56,177  0/ 

0014  P  A  I  A  I  TSRR/53  I  , I  023,636,  102, 0/ 

0  015  D  A I  A  SKUR/2  8  «  ,  -  I  0  78  ,1  285,0  94.-5  76/ 

0016  LU01=i 

001  7  NUf1  =  5 

0018  AC0R-AU/10  0 

0019  I  T  S  S  = I  TS/100 

0020  600  00  6  10  I  *  I  ,  NUN 

0  021  C -<  COVRt  I  >»  <24<<TC0RRC  I  >/IO  H*2*<  I  T$RR<  1  )/lO 0  )**2 

0  02  2  1  « ( SKUR<  I  )  ♦  I  0  )  •  *2  )/2 

0023  D 1 S  <  I  >=<  COV*COVR< I  >*ACOR«<  ACORR< 1  )/l 0  )  • I ISS* 1 1 SRR< I  )/ 100 

002  4  l*  (  SKWMC  >♦<  $KWR<  I  )♦  10  )  >-C 

0025  610  CONTINUE 

002  6  I C  L  A  S -  1 

002  7  0  !CP=D  l$<  1  > 

0028  DO  620  J*2,NUH 

0029  I  f ( D  I  S  <  J  >  CT  8  I  C  D  >  ICLAS*J 

0030  I  F  <  D  I  S ( J  )  GT  eiGO)  BICO^DIS(J) 

003  1  62  0  C  OH  1  i  HUE 

0032  GO  TO  ( 630, 640 .650 , 660, 670 >  ICLAS 

003  3  630  WR 1 1 E< LUO  T ,631) 

0034  GOTO  800 

0035  64  0  WRI I E( LOOT , 64  I  ) 

0036  GOTO  800 

003  7  650  WRI  TE<  LOOT, 651  > 

0038  GOTO  800 

0039  660  WR  1  TE<  LUO  1,661) 

004  0  G  U 1  0  800 

004  1  6  70  UR  I  TE<  LUOT  ,  67  1  ) 

0042  GGT0800 

004  3  63  1  F  OR  HA  T  <  *  C  I  T  Y*  ) 

004  4  64  l  F  OR  H  A  T  < *F  IELD*  ) 

^015  651  F  OR  H  A  I  < "FARNLAHD") 

004  6  66  1  F  ORNATE “FOREST  ") 

0047  671  FORMATE ‘WATER  *> 

004  8  800  IFELUOT  EO  6)  WR I T E ( L UO T  ,  8 50  ) 

0049  850  r ORNATE • I • > 

0050  RETURN 

0051  END 

0052  END! 
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